Equivalence between the mobility edge of electronic transport on disorder-less 
networks and the onset of chaos via intermittency in deterministic maps 
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We exhibit a remarkable equivalence between the dynamics of an intermittent nonlinear map and 
the electronic transport properties (obtained via the scattering matrix) of a crystal defined on a 
double Cayley tree. This strict analogy reveals in detail the nature of the mobility edge normally 
studied near (not at) the metal-insulator transition in electronic systems. We provide an analytical 
expression for the conductance as function of system size that at the transition obeys a (/-exponential 
form. This manifests as power-law decay or few and far between large spike oscillations according 
to different kinds of boundary conditions. 



PACS numbers: 05.45.Ac, 71.23.An, 71.30.Bs,47.52.+j 

Occasionally the detection of a deep running analogy 
between two apparently different physical problems al- 
lows for the determination of elusive quantities and un- 
derstanding of difficult issues. Here we present a re- 
lationship between intermittency and electronic trans- 
port. This development brings together fields of research 
in nonlinear dynamics and condensed matter physics. 
Specifically, the dynamics at the onset of chaos appears 
associated to the critical conductance at the mobility 
edge of regular self-similar networks [l[ . 

Recently, the dynamics at the transitions to chaos that 
occurs along the three known universal routes from reg- 
ular to irregular behavior (in low-dimensional nonlinear 
maps) has been analyzed with a good deal of detail JU- 
JU. This effort has helped establish the nature of the 
statistical-mechanical structure obeyed by the dynamics 
associated to nonmixing and nonergodic attractors J4j. 
On the other hand, there are known connections between 
nonlinear dynamical systems and electronic transport 
properties. For example, there are models for transport 
in incommensurate systems, where Schrodinger equations 
with quasiperiodic potentials [f| are equivalent to nonlin- 
ear maps with a quasiperiodic route to chaos, and where 
the divergence of the localization length translates into 
the vanishing of the ordinary Lyapunov coefficient Q . 

At the tangent bifurcation Q , the focal point of the in- 
termittency route to chaos, an uncommon but welcome 
simplicity has led to analytical results in closed form for 
the dynamics at vanishing Lyapunov exponent Q. Here 
we make full use of this circumstance showing that trans- 
port in a model network, a double Cayley tree, resolved 
by means of the scattering matrix, is given by the prop- 
erties of a one-dimensional nonlinear map. The model, in 
this study, does not contain disorder; nevertheless it dis- 
plays a transition between localized and extended states. 
The translation of the map dynamical features into elec- 
tronic transport terms provides not only the description 
of the two different conducting phases but, we believe, 
offers for the first time a rigorous account of the conduc- 
tance at the mobility edge. A type of localization length 



in the incipient insulator mirrors the departure from ex- 
ponential sensitivity to initial conditions at the transition 
to chaos. 

We recap briefly the usefulness of Cayley tree networks 
in the study of electronic transport properties in the pres- 
ence, and absence, of disorder. A single Cayley tree spans 
over a space of infinite dimensionality [7[ and transport 
on it exhibits a metal-insulator transition as a function of 
disorder Q. A scattering approach was applied in Ref. 
9] for off-diagonal disorder and shown that the metal- 
insulator transition occurs for connectivity K > 2 (K + 1 
is the coordination number). A single Cayley tree is a 
first approximation to an ordinary regular lattice but, 
as shown below, a double Cayley tree (two single Cayley 
trees joined conformally as in Fig. [T]) is a much better 
approximation (see also [l(|). In Ref. is shown that 
the conducting band of a disorder-less double Cayley tree 
contracts and disappears as K increases. In Ref. [12j the 
dynamic behaviour of a chain of scatterers was analyzed 
in the absence and presence of disorder; while the local- 
ization transition for different types of complex networks, 
including the double Cayley tree, was studied in Ref. [l3[ 
via spectral statistics. Finally, the double Cayley tree 
problem is relevant for transport in chaotic cavities with 
broken mirror symmetry [14| . Here we show that the 
electronic structure for K — 2 can be determined by re- 
ducing the scattering matrix to a nonlinear map. This 
development facilitates the band description of the con- 
ductance as a function of energy including the location 
of the mobility edge. 

Here we consider electronic transport in the double 
Cayley tree (see Fig. [T]). We refer only to the ordered, 
crystal-like, system and reduce its associated scattering 
matrix to a nonlinear map. The number of times the trees 
are ramified, starting from perfect join, is the generation 
n that quantifies the size of the system. Also, for brevity, 
we will fix the tree connectivity to K = 2 where one lead, 
we call it the incoming lead, is divided into two leads at 
a given node. The leads are assumed to be equivalent to 
one dimensional perfect wires with length equal to the 
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FIG. 1: A double Cayley tree of conectivity K = 2 and lattice 
constan a. Each bond is a perfect one dimensional conductor. 



lattice constant a and are independent of n. Hence, each 
node is described by the same 3x3 scattering matrix for 
which we assume the model 
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where e, a real number in the interval < e < i is 
the transmission probability (or coupling) from the in- 
coming lead to the others, and viceversa. The reflec- 
tion amplitude to the incoming lead is — (a + f3), with 
a = -(1 - VT~^2i)/2 and /3 = (1 + v / l _ ^2e)/2. When 
incidence is only on one of the other two leads a is the 
reflection amplitude to the same lead and (3 the trans- 
mission amplitude to the other lead. 

The scattering matrix of the system is 2 x 2 and sat- 
isfies a recursive relation. If we denote by S n the scat- 
tering matrix at generation n, the combination rule for 
scattering matrices allows S n to be written in terms of 
the scattering matrix at a previous generation n — 1, 
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with 1 the 2x2 identity matrix. The scattering matrix at 
a generation n can be obtained iteratively starting from 
that for the perfect union in the middle of our double 
Cayley tree: So — <J X , where a x is a Pauli matrix. 

Firstly, it can be seen that S n is a unitary matrix, 
which is the condition of flux conservation. Then, time 
reversal invariance restricts S n to be a symmetric matrix. 
Finally, the additional lattice spatial reflection symmetry 
implies that S n has the form |14| 



S n — 



(3) 



where r n and t n are the reflection and transmission am- 
plitudes. With this structure S n is diagonalized by a 
tt/ 4- rot at ion that is 
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Here, n and 9' n are the eigenphases that satisfy e l9n = 
Tn+Ui and e % ® n = r n — t n - In terms of the eigenphases the 



(i.e. 



transmission amplitude is given by t r 
Morever, the dimcnsionlcss conductance (i.e. in units of 
2e 2 /h) depends on the eigenphases through Landauer's 
formula as [H,[i3, 



fjn = \t n f 



(5) 



Therefore, the analysis of the eigenphases is of crucial 
importance as they determine every transport property. 

Central to our discussion is the fact that the recursive 
relation @ can be written in the diagonal form ((J) and 
this implies the existence of a one-dimensional nonlin- 
ear map for the phase 6 n . The map 6 n +i = f{9 n ) can 
actually be obtained in the following closed form 
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where the dependence on e and ka comes out clearly. In 
what follows everything said about n is valid for 0' n as 
well. Perfect union at n — means 9o = and 9' Q = tt. 

For a given value of e the map ([6]) is periodic in ka 
(the parameter related to the energy) with period n. In 
the range < ka < tt the attractor diagram presents a 
chaotic region between two windows of period 1 (see inset 
of Fig. d{b)), separated by bifurcation points at 
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arccos \/2e, (7) 



such that the chaotic region of the map takes place in 
the interval k c a < ka < k c ia, while windows of period 
1 in ka < k c a and ka > k c *a. As we see below Eq. ([7]) 
gives the locations of the mobility edge as a function of 
the transmission probability e. Fixed-point solutions #oo 
for 6„ , n — > oo , are 



9_ for < ka < k c a 
1 for k c a < ka < k c >a 
9 + for k c ia < ka < tt 

where 8 = —ka + 3tt/2, and 9± is given by 
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tan 9 4 



sin ka (cos ka ± Vcos 2 ka — 2e) 



1 — cos ka (cos ka ± a/ cos 2 ka — 2e) 



(9) 



These fixed-point solutions indicate that for large n, 9 n 
reaches the values 9± in the windows of single period, 
while in the chaotic region 9 n fluctuates according to an 
invariant density with maximum at 9. 

At the bifurcation points k c a and k c >a, the fixed-point 
phase 9oo takes the critical values 9 C = — arccos \/2e + 
3tt/2 and 9 C > = arccos \J~2t + tt/2, respectively, for 9oo 
between and 2-7T. In Fig. [21a) we plot the map 
Eq. ^ for e = 1/4 for the three parameter values 
ka = 7r/4, 27r/4, 37t/4. It is evident that the transitions 
from chaotic behavior to the windows of period one are 
tangent bifurcations as the map is tangent to the identity 
line at the critical values k c a — tt /A and k c >a — 37r/4, 




FIG. 2: f(6 n ) of Eq. © for e = §: (a) fca = tt/4 (dashed), 
27r/4 (long dashed), 37r/4 (continuous), and (b) ka = 2.6 
(long dashed), 37r/4 (continuous), 2.1 (dashed). The dotted 
lines correspond to the identity. Inset: Periodic and chaotic 
attractors. 
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FIG. 3: Conductance as a function of generation for e = 4: 
Am = 0.5 (left panel) and ka — ir/i (right panel). Continu- 
ous lines represent g n obtained directly from the map through 
Eq. ([5]) while dotted lines correspond to g n = exp(2Ain) and 
g n = exp 3 y 2 (2A3/2n) for left panel and right panel, respec- 
tively. The two curves in (b) differ because of a proportion- 
ality factor, See Eq. (14). 



where 9 C — 5ir /4 and 9 C > = iir/A. Fig. [^b) shows three 
cases: secant (ka — 2.6) and tangent (ka = 37r/4) period 
one solutions, and a bottleneck (ka = 2.1) that gives rise 
to intermitency as a precursor to periodic behavior. 

Information about transport can be obtained from the 
sensitivity to initial conditions that characterizes the dy- 
namics of the nonlinear map. For finite n it is given by 



d9 r , 
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exp(Ai(n) n), 
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where 9q is an initial condition and the exponential law 
after the 2nd identity defines the finite n Lyapunov ex- 
ponent Ai(n). For n large Ai(n) becomes the Lyapunov 
exponent Ai, a number independent of 9q that according 
to its sign characterizes periodic and chaotic attractors. 
At the tangent bifurcation Ai = and the sensiti vity 
adopts instead a g-exponential form (see below) 0, lle| . 
FromEqs. ©, (TQJ and Eq. © (and t n = ±(e lS "-e<)) 
we obtain the recursion formula 

g n = ffn-i exp(Ai(n) n) exp(A' 1 (n) n), (11) 



Ax(n) = In 



1 - e - Vl - 2ecos(6»„ + 2ka) 
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where A.[(n) is given by Eq. (fT2")) with 9 n replaced by 9' n . 
We note that Ai and A[, and hence g n , do not depend on 
the initial conditions 9q and 9' . The Lyapunov exponent 
Ai is given by Eq. (T2"j) with 6 n replaced by 9^. 

The dynamical properties of the map © translate into 
the following network properties: i) In relation to the at- 
tractors of period one that take place along < ka < k c a 
and k c >a < ka < ir, wc corroborate from Eqs. ([5]) and 
(fT2")l that Ai is negative and therefore the conductance g n 
decays exponentially with system size n, g n — exp(2Ain), 
implying localization, the localization length being £i = 
a/|Ai|. In the left panel of Fig. Owe see a clear exponen- 
tial decay of g n as a function of n at ka = 0.5, where we 



compare g n computed directly from Eq. ([5]) with that ob- 
tained from Ai. We notice that the conductance for the 
localized states of an ordered system displays the same 
behavior as that in the insulating regime of a disordered 
wire in quasi-one-dimensional configuration |19| . ii) With 
respect to the chaotic attractors that occur in the inter- 
val k c a < ka < k c /a we observe that Ai becomes positive 
and the recursion relation Eq. (fTTj) does not let g n decay 
but makes it oscillate with n (not shown here) indicat- 
ing that conduction takes place. In our model g n does 
not scale with system size as in the metallic regime of 
quasi-one-dimensional disordered wire where Ohm's law 
is satisfied 19] . In the parameter region where the map is 
incipiently chaotic, say ka > k c a, the network grows with 
n with an insulator character, but interrupted for other 
intermediate values of n by conducting crystals. In the 
map dynamics these are the laminar episodes separated 
by chaotic bursts in intermittent trajectories. 

The most distinct outcome of our treatment is the de- 
scription obtained of the mobility edge from the dynam- 
ics at the critical attractors located at ka = k c /a and 
ka = k c a. There Ai = and according to Eq. (fTTj) 
not much can be said about the size dependence of the 
conductance when n> 1. However we can use to our ad- 
vantage the known properties of the anomalous dynamics 
occurring at these attractors once they are identified as 
tangent bifurcations [see Fig. f^b) for ka = 3tt/4\. At 
a tangent bifurcation of general nonlinearity z > 1 the 
sensitivity obeys a ^-exponential law for large n 0, 

£„ = H n> i = exp 9 (A«jn) = [l-(q- l)\ q n] T ~, (13) 

where X q is a ^-generalized Lyapunov coefficient given by 
X q = ^fzu, q = 2 — 1/z, where u is the leading term 
of the expansion up to order z of 9 n close to 9 C (or 
9c); i.e. 8 n - 9 C = (0 n _r - 9 C ) + u - 9 C \ Z + ■■■. 
The minus and plus signs in Eq. (|13p and in X q corre- 
spond to trajectories at the left and right, respectively, 
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of the point of tangency 9 C . Eq. (fl~3f implies power- 
law decay of £ n with n when 9 n — 9 C < and faster 
than exponential growth when 8 n — 9 C > 0. (We recall 
that any choice of So other that the Pauli matrix trans- 
lates into another initial condition for the map). By 
making the expansion around 9 C (or 6 c i) for our map 
(|o| we find (as evidently anticipated) z — 2 implying 
q = 3/2, u = \J (1 — 2e)/2e, and the ^-generalized Lya- 
punov exponent is A3/2 = —2^/(1 — 2e)/2e. Following 
the same steps that lead to Eq. (fTT]) . the recursion re- 
lation for g n at each bifurcation point takes the form 
9n = 5™-iexp 3/2 (A 3/2 (n)n)exp 3/2 (A 3/2 (n)n), so that 
when 9 n — 9 C < (see Eq. (flB")) 1 ), g„ = exp 3 / 2 (2A 3 / 2 n), or 

9n oc ^1 - ^A 3/2 n^ . (14) 

In the right panel of Fig. [3] we compare the results from 
Eq. ([5]) (continuous line) and Eq. (fl4f (dashed line). It 
is clear that g n decays as a power law (with quartic ex- 
ponent) rather than the exponential in the insulating 
phase. We enphasize that a localization length given by 
C3/2 = a /^3/2 can still be defined at the mobility edge. 
To our knowledge this property has not been reported 
before. When 9 n - 9 C > 0, A 3/2 = 2^(1 - 2e)/2e > 
and the recursion relation for g n describes, as the re- 
sult of the diverging duration of the laminar episodes of 
intermittency, large n intervals of vanishing g n between 
increasingly large spike oscillations. 

In summary, we can draw significant conclusions about 
electronic transport from our study. These arise natu- 
rally when considering the dynamical properties of the 
equivalent nonlinear map near or at the intermittency 
transition to chaos. Since iteration time in the map trans- 
lates into the generation n of the network, time evolu- 
tion means growth of system size, reaching the thermo- 
dynamic limit (and true self similarity) when n — > 00. 
In that limit, windows of period one separated by a 



chaotic band correspond, respectively, to localized and 
extended electronic states. Further, in the referred pa- 
rameter (feet, e) regions the conductance g n of the model 
crystal shows either an exponential decay with system 
size, with localization length given by £1 (&s in the case 
of a quasi-one-dimensional disordered wire 19]), or an os- 
cillating property signalling conducting states. The pair 
of tangent bifurcation points of the map correspond to 
the band or mobility edges that separate conductor from 
insulator behavior. At these bifurcations the sensitivity 
to initial conditions exhibits either power-law decay 
(when #0 < $c' or 9q > 9 C ) or faster than exponential 
increase (9 c i < 9q < 9 C ) and consequently the conduc- 
tance inherits comparable decay or variability with sys- 
tem size n. Notably, as we have seen we can still de- 
fine a localization length, the g-generalized localization 
length £ 9 with a fixed value of q = 3/2. This expression 
is universal, i.e. it is satisfied by all maps that in the 
neighborhood of the point of tangency have quadratic 
term, i.e. z — 2 [2|. This quantity can be obtained di- 
rectly by evaluation of g n (in Eq. ([5])) when n — ► 00. 
At the mobility edge C-f 1 = — linin^oo n^ 1 lng n van- 
ishes because In g n no longer decreases linearly with n, 
as it is the case in the insulating phase. However, use of 
Q -1 = — linin^oo n^ 1 ln 9 g n leads to a finite number for 
one particular value of q, q = 3/2, when the degree of de- 
formation q in the y-logarithm restores linear behavior. 
Otherwise £ 9 vanishes or diverges. 

In spite of the unusual features of the double Cayley 
tree transport model, the complete set of exact solutions 
derived from it provides a comprehensive picture about 
non-exponential behavior of central quantities like the 
conductance at the transition between the insulator and 
conductor regimes. 
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